Decoherence due to elastic Rayleigh scattering 
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We present theoretical and experimental studies of the decoherence of hyperfine ground-state 
superpositions due to elastic Rayleigh scattering of light off-resonant with higher lying excited 
states. We demonstrate that under appropriate conditions, elastic Rayleigh scattering can be the 
dominant source of decoherence, contrary to previous discussions in the literature. We show that 
the elastic-scattering decoherence rate of a two-level system is given by the square of the difference 
between the elastic-scattering amplitudes for the two levels, and that for certain detunings of the 
light, the amplitudes can interfere constructively even when the elastic scattering rates from the 
two levels are equal. We confirm this prediction through calculations and measurements of the total 
decoherence rate for a superposition of the valence electron spin levels in the ground state of Be + 
in a 4.5 T magnetic field. 



Off-resonant light scattering (spontaneous emission) is 
an important source of decoherence in many coherent- 
control experiments with atoms and molecules. Exam- 
ples include the use of optical-dipole forces for gates in 
quantum computing |l[ , the generation of spin squeezed 
states through laser-mediated interactions [2|-|6[ , and the 
trappingand manipulation of neutral atoms in optical 
lattices 0, Q • These experiments frequently involve su- 
perpositions of two-level atomic systems (qubits) and use 
laser beams off-resonant with higher lying excited states 
to control and measure the atomic states. 

In general, decoherence of an atomic superposition 
state due to off-resonant light scattering occurs if the 
scattered photon carries information about the qubit 
state. During Raman scattering the initial and final 
qubit states differ. The state of the scattered photon 
is entangled with the atomic state, providing "welcher- 
weg" (which-way) information and leading to decoher- 
ence [9J, [lCf . By contrast the role of elastic Rayleigh scat- 
tering for decoherence is not as clear. Two very different 
regimes have been discussed and are supported by ex- 
periment. On the one hand it has been found that in 
some experiments Rayleigh scattering gives rise to neg- 
ligible decoherence provided that the elastic scattering 
rates from both qubit levels are approximately equal [lOj ■ 
On the other hand, decoherence due to Rayleigh scatter- 
ing of photons on a cycling transition is used for strong 
projective state measurement 11|. 

In this letter we develop a microscopic theory for the 
decoherence of a qubit due to elastic Rayleigh scat- 
tering that gives a unified treatment of these different 
regimes. Our key finding is that the decoherence induced 
by Rayleigh scattering is proportional to the square of the 



difference of the probability amplitudes for elastic scat- 
tering from the two levels [Eq. |0]. When the two ampli- 
tudes are approximately equal the resulting decoherence 
rate can be small (first case above) and when one ampli- 
tude dominates the other, Rayleigh decoherence can be 
large (second case above). However, even when the elas- 
tic scatter rates are approximately equal, the scattering 
amplitudes can have opposite sign - a situation appli- 
cable to many coherent control experiments [2|-|6| - and 
Rayleigh decoherence can dominate Raman decoherence. 

In the rest of this paper we derive the microscopic 
model that quantitatively describes the decoherence due 
to Rayleigh scattering and verify our predictions against 
an experiment with 9 Be in a 4.5 T magnetic field. For 
concreteness we discuss the problem in terms of the rele- 
vant energy level structure of 9 Be + shown in Fig. [TJ The 
two ground state sub- levels \u) = \ + i) and \d) = \ — |) 
corresponding to the valence electron spin states parallel 
and anti-parallel to an applied magnetic field are split 
by il z = 2ttx 124.1 GHz. Also shown are the first ex- 
cited P state levels \J, Mj), J = |,|. We consider a 
non-resonant light field (X)a=-i &a£a)£o cos(ko -x — Wot), 
where b\ is the amplitude of the corresponding polar- 
ization component i\ of the light, and derive a master 
equation for the time evolution of the | ± ~) ground state 
sub- levels in the presence of the off- resonant light. 

An effective Hamiltonian coupling of the two qubit lev- 
els via spontaneous scattering of the driving field into the 
vacuum modes can be derived by adiabatic elimination 
of the excited states shown in Fig. Q] and is given by 



H sp — Hdd 



H u 



H, 



(It I. 



H U a 



(i) 



'Electronic address: huys@csir.co.za 
^Electronic address: john.bollingerSnist.gov 



Here the different terms describe processes which either 
do, or do not, flip the atomic spin upon scattering of a 
photon, as indicated by the subscripts. The Raman terms 
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FIG. 1: Relevant energy level structure of 9 Be + at 4.5 T. We 
only show the mi — + 1 levels which are prepared experimen- 
tally though optical pumping. The off-resonant laser is shown 
with a frequency uj appropriate for the experimental results 
presented in Fig. [3] 



that lead to spin- flips are (in the Schrodinger picture), 

H du = hJ2 (h d A x a-a kx E-(Ak) + h d ,Ha + al x E + (Ak^ 

JkA 

H ud = hY / (h u jl x a-ai x E + (A]s.) + ^* A «T + a kA £-(Ak)) 



where a fcA (cikA) is the raising (lowering) operator for the 
k, A mode of the electromagnetic field, a + — (a x + ia y )/2 
and a~ = (a x —ia y )/2 where a x , a y , and a z are the Pauli 
matrices for the two- level ground state, and E + (Ak) = 
J5T*(Ak) = (£ /2)e t(ko - k) ' x - l " Jot . The coupling constants 
are 
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where i, j E {-±, |}, 5u\ = £lj.\+i - Oj - w is the 
detuning of luq from the |i) — > |J, Mj = A + i) tran- 
sition, at energies Mli and hVlj\ + i respectively, and 

3k = \/hbJk/(2eoV) with V the quantization volume. Sim- 
ilarly the Rayleigh terms are 



H dd = hJ2 [hf kx a kx E- (Ak) + h d r i\at x E+(Ak)j <h 

JkX 

H uu = hJ2 (h U jlA x E-(Ak) + hT k lal x E+(Ak)) a 2 , 



where we've defined the operators &i = \{I — <r z ) — 
\d)(d\, d"2 = \{I + u z ) = \u)(u\. With these definitions 
the Rayleigh terms can be combined H e i = Hdd + H uu 
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The first term in Eq. ([2]) is proportional to the identity 
for the qubit and consequently does not lead to decoher- 
ence. The second term proportional to a z leads to pure 
dephasing which, as we will see, can sometimes be the 
dominant decohering process. 

The time evolution of the reduced density matrix, 
Ps{t'), of the qubit is well described by the master equa- 
tion [TJ] (in the interaction picture) 

^1^ = -± £ dt'Trk{[H sp (t),lH sp (t'),p S (t')p*]}}, 

(3) 
in which ps(t') is the density matrix of the qubit before 
tracing over all vacuum modes represented by density 
matrix p k as indicated by Tr k . 

In evaluating Eq. Q we find_Raman scattering terms 
given in the Lindblad form by [19f 
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Tij is the rate for an ion initially in state \i) to scatter 
a photon and end up in state \j) and is given by the 
Kramers-Heisenberg formula 
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Here VLr = p£o/(2h), 7 is the spontaneous decay rate of 
the excited states, and the amplitudes A z y^ for scattering 
through intermediate level | J, A + i) are defined as 

. _ fex(j|d • e£ +(w) | J, A + i)(J, A + »|d • ex\i) 

A , \ — ; (5) 
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where p = (|, ||d • £1)5) is the largest dipole matrix 
element. We define the total Raman scattering rate 

t Ram t du 1 t u d~ 

The contributions to Eq. ([3]) from Rayleigh scattering 
are 

Cdd,uuPs(t) = — |- \ps(t) - cf z ps{t)a z ] (6) 



where 
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Finally the full time evolution of the reduced density 



matrix is 
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Eq. ([7]) is the main theoretical result of this letter. 
General master equation treatments of decoherence due 
to elastic scattering have appeared in the literature (13J . 
However, to our knowledge this is the first calculation 
producing a simple formula which demonstrates the cen- 
tral role of elastic scattering amplitudes in spontaneous- 
emission- induced decoherence of hyperfine-state superpo- 
sitions. 

Since the detunings in the two amplitudes in Eq. ([7]) 
can have opposite sign, the different contibutions can add 
constructively, leading to a large Rayleigh scattering de- 
coherence. Physically the terms proportional to a z in 
Eq. @ suggest that this decoherence is due to dephasing. 
A photon elastically scattered into mode k, A produces a 
shift in the phase of a superposition of the |u), \d) levels 
proportional to the difference in the amplitudes for scat- 
tering into that mode. Since photon scattering occurs 
at random times and into random directions, random 
phase shifts are produced (due to the _E + (Ak) coefficient), 
and consequently phase information of the superposition 
state is lost. 

We realize elastic scattering amplitudes with opposite 
sign in a Penning ion trap operating in the strong-held 
Zeeman regime. These traps have found increasing im- 
portance in quantum information studies 14H16J] , under- 
scoring the need for quantitative analysis of decoherence 
processes in these systems. Details of our basic Pen- 
ning trap set-up can be found in [17|. We use two 313 
nm laser beams (see Fig.QJ to Doppler laser cool (T~0.5 
mK) and optically pump a planar crystal (diameter —300 
/im) of -100 9 Be+ ions into the |u) = | + 1/2} state. A 
phase- locked Gunn diode oscillator near 124 GHz pro- 
vides microwave radiation used to perform high-hdelity 
(>99 %) qubit rotations. At the end of an experiment, 
the population in the \u) state is measured by detecting 
the resonantly scattered ion fluorescence induced by the 
Doppler cooling laser (tuned to the \u) — > |3/2, 3/2) cy- 
cling transition), and normalizing this fluorescence level 
to that obtained with all the ions in \u). 

We introduce an off- resonant, linearly polarized laser 
beam and measure the total decoherence rate | (Tji am + 
T e i) from a decrease in length of the composite Bloch vec- 
tor of the ions using a two-7r-pulse spin echo sequence like 
that shown in Fig. [2£a). We adjust the relative phases 
of the 7r/2 and ir pulses so that in the absence of any 
coherence loss all ions are rotated to \d) at the end of 
the spin-echo sequence. The measured population in \u) 
therefore provides a measure of the total decoherence. 
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FIG. 2: (a-c) Schematic pulse sequences used to measure (a) 
total decoherence and (b), (c) Raman scattering rates. Ra- 
man scatter rates Fd u and F u d (used to calibrate laser inten- 
sity, see Fig. [3} obtained by initializing the ions in \d) or \u) 
respectively and measuring the |m) state population vs off- 
resonant laser interaction time, (d) Example measurements 
of the total decoherence and Raman scatter rates. Total deco- 
herence (solid markers), left axis; Raman scatter rates (open 
markers), right axis. 



From Eq. (|8]) we calculate that at the end of the spin 
echo sequence p uu = -|(1 — e ~( rR '"™ +rc! ) r / 2 ), where r is 
the total laser interaction time of the spin echo sequence. 
Figurej^d) shows measurements of the \u) state popula- 
tion as a function of the spin echo laser interaction time. 
We determine \(TR am + r e ;) from twice the fitted slope 
at short spin-echo sequence times. 

Other sources of decoherence can be neglected. With- 
out off-resonant light we measure a — 1 % loss in co- 
herence due to magnetic held fluctuations for up to 6 ms 
spin echo times, much longer than the spin echo sequence 
times used to determine decoherence due to off-resonant 
light scattering. Further, the ac-stark shift from the off- 
resonant laser is nulled to better than 1 kHz by rotating 
the polarization. Doing so minimizes decoherence due to 
laser-power fluctuations during the spin-echo probe se- 
quence. Electric field homogeneity due to the laser beam 
waist is greater than 95 % over the ion planar array. 

Figure 02a) presents two sets of measurements of the 
total decoherence rate as a function of the frequency of 
the off-resonant light, and compares these measurements 
with different theoretical predictions. For these mea- 
surements the intermediate state detunings and result- 
ing d — >■ d and u — > u scattering amplitudes in Eq. ([7]) 
have opposite sign, producing a large elastic decoherence 



rate. This condition is in part achieved because the light 
detunings to the intermediate -P3/2 levels are comparable 
to the large (124 GHz) qubit splitting. 

The two measurement sets use different ion crystals 
and different techniques for calibrating the off-resonant 
light intensity (see Figs[3Jb),(c)). Good agreement is ob- 
tained with the full theory presented here which man- 
ifests a Rayleigh decoherence rate determined by the 
square of the difference between the \u) and \d) elastic 
scattering amplitudes (Eq. ([?[) )■ This contrasts with the 
idea conveyed in previous literature that decoherence due 
to elastic scattering is produced by a difference in elas- 
tic scattering rates. For example Ref. [18J] (Sec. III.B) 
discusses the which-way information acquired due to a 
difference in Rayleigh scattering rates and estimates a 
decoherence rate ^T e i diff — \ "" , * d \ due to these rate 

differences. Figure EJa) shows theory based on Y e i t< iiff- 
As the laser detuning approaches resonant transitions 
at —79.4 GHz and —37.7 GHz, reasonable agreement 
is obtained with this theory. However, at a detuning 
of ~ —56 GHz, Y uu « Tdd and we measure a decoher- 
ence rate that is five times larger than that predicted by 

These measurements can be contrasted with previous 
work [lCf using trapped 9 Be + ions at low magnetic field 
in which Rayleigh-scattering-induced decoherence was 
shown to be negligible. The results of that experiment 
benefited from some fortuitous conditions: the detuning 
A of the light from any atomic excited state was large 
compared to the qubit transition frequency, and the qubit 
states were clock states (states whose energy difference 
does not depend to first order on changes in the magnetic 
field). These conditions imply a nearly complete cancella- 
tion of the elastic scattering amplitudes in Eq. (JTJ) result- 
ing in a weak Rayleigh decoherence with a dependence on 
detuning ex 1/A 4 . For general two-level systems the can- 
cellation will not be as complete, resulting in a Rayleigh 
decoherence rate oc 1/A 2 . The prescription discussed in 
this manuscript now enables an accurate calculation of 
Rayleigh decoherence for these low-field trapped ion as 
well as other coherent control experiments. 
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